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Sensitivity  of  Satellite  Multichannel  Sea  Surface  Temperature 
Retrievals  to  the  Air-Sea  Temperature  Difference 

Douglas  A.  May  and  Ronald  J.  Holyer 


Naval  Research  Laboratory,  Remote  Sensing  Applications  Branch,  Stennis  Space  Center,  Mississippi 

The  advanced  very  high  resolution  radiometer  multichannel  sea  surface  temperature  (MCSST) 
retrieval  technique  provides  global  algorithm  accuracy  statistics  generally  showing  a  bias  of  less  than 
0. 1°C  and  an  rms  error  of  less  than  0.7°C  when  compared  to  colocated  drifting  buoy  in  situ  data  in  the 
absence  of  aerosols.  This  remaining  error  is  not  always  random  but  is  shown  to  be  correlated  to  the 
air-sea  temperature  difference.  The  MCSST  technique  is  modeled  and  then  compared  to  in  situ  data  to 
show  this  dependency.  Atmospheric  radiative  transfer  calculations  are  used  to  provide  a  simulation  of 
satellite  retrieval  sensitivity  to  air-sea  temperature  differences.  Buoy  sea  surface  temperature  (SST) 
and  air  temperature  observations  are  then  presented  as  experimental  verification  of  the  simulation 
results.  Retrieval  errors  depend  both  on  the  mean  air-sea  temperature  difference  conditions  present  in 
the  data  set  used  to  empirically  derive  the  algorithm  and  on  the  changes  in  air-sea  temperature 
difference  conditions  relative  to  the  derivation  data  set  mean  conditions.  Retrieval  error  is  found  to 
respond  linearly  with  air-sea  temperature  difference  changes.  MCSST  retrieval  errors  of  1.0°C  can 
occur  for  air-sea  temperature  difference  changes  of  7°-10°C  from  mean  conditions  when  the 
dual-window  (channels  3  and  4)  or  triple-window  (channels  3,  4,  and  5)  algorithms  are  used.  The 
split-window  (channels  4  and  5)  MCSST  algorithm  is  shown  to  be  less  sensitive  to  air-sea  temperature 
differences.  Cross-product  SST  (CPSST)  and  nonlinear  SST  (NLSST)  algorithms  are  also  examined. 
These  algorithms  generate  results  similar  to  the  MCSST  algorithm  for  the  dual-  and  triple-window 
equations.  However,  the  CPSST  and  NLSST  split-window  algorithms  demonstrate  greater  sensitivity 
to  air-sea  temperature  difference  changes  than  do  the  MCSST  split-window  algorithm.  Retrieval  errors 
of  1°C  can  occur  for  air-sea  temperature  difference  changes  of  W-n'C  from  mean  conditions.  Users 
of  satellite  SST  retrievals  in  regions  that  experience  large  fluctuations  in  air-sea  temperature  difference 
should  be  aware  of  this  possible  error  source. 
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Introduction 

The  generation  of  absolute  sea  surface  temperature  (SST) 
from  satellite  infrared  (IR)  brightness  temperature  measure¬ 
ments  has  been  studied  often  [Andtrtg  and  Kauth,  1970; 
Maul  and  Sidran,  1972;  McMillin,  1975;  Bernstein,  1982; 
Llewellyn- Jones  et  al.,  1984;  McMillin  and  Crosby,  1984; 
Walton,  1988;  Barton  and  Cechet,  1989;  Barton  et  al.,  1989]. 
Because  the  surface-emitted  radiation  is  partially  absorbed 
and  reradiated  by  the  atmosphere,  an  atmospheric  correc¬ 
tion  must  be  made  to  the  satellite-measured  radiances  before 
the  (SST)  can  be  determined.  Typically,  two  or  more  IR 
measurements  at  differing  wavelengths  or  viewing  angles  are 
used  to  determine  this  atmospheric  correction  [McMillin  and 
Crosby,  1984]. 

Global  multichannel  sea  surface  temperature  (MCSST) 
retrievals  from  satellites  have  been  operationally  produced 
from  National  Oceanic  and  Atmospheric  Administration 
(NOAA)  advanced  very  high  resolution  radiometer 
(AVHRR)  data  for  several  years  at  the  National  Environ¬ 
mental  Satellite  Data  and  Information  Service  (NESDIS) 
[e.g.,  McClain,  1989].  MCSSTs  are  retrieved  globally  each 
day  on  an  orbit-by -orbit  basis  from  global  area  coverage 
(GAC)  data  at  a  spatial  resolution  of  8  km  (2  by  2  GAC 
arrays).  Details  about  operational  MCSST  processing  and 
cloud-screening  algorithms  are  given  by  McClain  et  al. 
(1985). 

Operational  retrieval  algorithms  are  empirically  derived 
by  regressing  satellite  retrievals  to  global  drifting  buoy  SST 
data  typically  within  4  hours  and  25  km  [Strong  and  Mc- 
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Clain,  1984].  This  is  routinely  accomplished  soon  after 
satellite  launch  and  rarely  requires  change  until  the  next 
satellite  launch.  Through  the  use  of  cloud-screening  algo¬ 
rithms  and  the  increased  quality  and  quantity  of  in  situ  data 
the  accuracy  of  satellite  data  has  improved  to  the  point  that 
present  monthly  global  comparisons  with  drifting  buoys 
generally  show  a  bias  of  less  than  O.TC  and  an  rms  error  of 
less  than  0.7°C  [McClain,  1989]. 

The  use  of  the  MCSST  technique  when  large  air-sea 
temperature  differences  exist  can  lead  to  reduced  accuracy 
[McMillin  and  Crosby,  1984].  Air  and  sea  surface  tempera¬ 
ture  difference  is  often  less  than  2'’C  [Isemer  and  Hasse, 
1985];  however,  coastal  locations  can  experience  greater 
extremes.  It  is  important  to  investigate  the  magnitude  of 
global  retrieval  algorithm  errors  under  such  anomalous  at¬ 
mospheric  conditions.  Accurate  knowledge  of  the  coastal 
water  environment  is  valuable  to  both  military  and  commer¬ 
cial  activities.  Since  coastal  locations  experience  frequent 
air-sea  temperature  extremes  not  often  observed  in  the  deep 
ocean,  identifying  satellite  retrieval  errors  under  such  con¬ 
ditions  will  enable  us  to  develop  improved  retrieval  tech¬ 
niques. 

Studies  using  simulated  retrieval  algorithms  have  pre¬ 
dicted  satellite  SST  errors  for  regionally  optimized  algo¬ 
rithms  when  used  under  anomalous  atmospheric  conditions 
[e.g.,  Minnett,  1986;  Minnett,  1990].  These  errors  were 
predicted  to  range  from  a  few  tenths  of  a  degree  to  more  than 
rc,  depending  on  the  algorithm  employed.  However,  such 
predictions  have  not  been  verified  with  in  situ  data. 

This  paper  investigates  the  sensitivity  of  multichannel  SST 
retrieval  techniques  to  air-sea  temperature  difference,  here¬ 
inafter  referred  to  as  ASTD.  A  simple  model  is  presented 
that  predicts  a  dependence  of  retrieval  errors  on  the  ASTD. 


CO 

tSD 

% 


12,567 


12,568 


.  May  and  Holyer:  Air-Sea  Temperature  Difference' 


The  LOWTRAN  7  [Kneizys  et  al.,  1988]  atmospheric  radia¬ 
tive  transfer  model  is  then  used  to  provide  a  simulated  error 
prediction  of  the  multichannel  technique  for  significant 
ASTD,  This  predicted  error  is  then  confirmed  utilizing 
coincident  satellite  SST  retrievals,  buoy  SST,  and  buoy  air 
temperature  observations.  In  addition  to  the  MCSST  algo¬ 
rithm  the  more  recent  cross-product  SST  (CPSST)  and 
nonlinear  SST  (NLSST)  algorithms  are  also  evaluated.  The 
implications  of  these  results  and  the  identification  of  condi¬ 
tions  under  which  satellite  SST  data  are  likely  to  be  less 
accurate  are  then  discussed. 

Theory 

If  one  assumes  that  the  atmosphere  is  a  horizontally  and 
vertically  homogeneous  absorbing  but  nonscattering  me¬ 
dium  and  that  the  emissivity  of  the  sea  surface  is  at  unity 
(implying  that  surface  reflectivity  is  zero),  the  following 
simple  radiative  transfer  equation  results: 


Fig.  I.  The  percent  of  transmittance  of  radiant  energy  as  a 
function  of  wavelength  for  a  wide  spectral  range  for  the  U.S. 
Standard  Atmosphere  model,  including  the  three  AVHRR  infrared 
channels  {LtiVioleUe,  1988]. 
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where  is  the  energy  received  at  the  satellite,  B(A,  T)  is 
Planck’s  function  for  energy  at  wavelength  A  emitted  by  a 
blackbody  at  temperature  T,  t(A)  is  the  transmittance  of  the 
atmosphere  at  wavelength  A,  is  the  temperature  of  the  sea 
surface,  is  the  mean  temperature  of  the  atmosphere,  and 
Al  and  A2  define  the  sensor  spectral  band  pass. 

Equation  (1)  describes  radiative  transfer  in  this  simple 
model  in  terms  of  radiant  energy.  The  radiative  transfer 
equation  expressed  in  terms  of  equivalent  blackbody  tem¬ 
peratures  results  in 


=  tT,  +  (1  -  t)T„  (2) 


where  ris  the  “eflfective”  value  of  7(A)  in  the  spectral  band 
between  A|  and  A2.  The  temperature  observed  by  the  satel¬ 
lite  sensor,  T^i,  is  equal  to  only  in  the  hypothetical  case 
where  t  =  1 . 

Equation  (2)  cannot  be  solved  directly  for  T,,  because  it 
contains  two  unknowns,  T,  and  t.  However,  two  simulta¬ 
neous  equations  can  be  obtained  by  sensing  the  ocean  in  two 
different  spectral  bands,  hence  the  multichannel  technique. 
A  complete  discussion  of  the  multichannel  solution  theory 
can  be  found  in  McMillin  [1975].  This  relationship,  ex¬ 
pressed  in  terms  of  channel  brightness  temperature  and 
atmospheric  transmittance,  results  in 
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where  T]  and  T-i  represent  the  IR  channel  brightness  tem¬ 
peratures,  respectively,  and  t,  and  T2  represent  the  atmo¬ 
spheric  transmittance  for  each  channel  bandwidth,  respec¬ 
tively. 

Figure  1  shows  an  atmospheric  transmittance  curve  appli¬ 
cable  to  the  LOWTRAN  7  U.S.  Standard  Atmosphere 
model.  The  spectral  bands  for  AVHRR  channels  3, 4,  and  5 
are  also  denoted  in  Figure  1.  Note  that  a  typical  transmit¬ 
tance  value  for  channel  3  is  approximately  0.9,  while  a 
typical  transmittance  value  for  channel  4  is  about  0.7.  If 
these  values  are  inserted  in  (3)  as  T|  and  t2,  the  following 
MCSST  algorithm  is  obtained: 


where  T(3.  4)  represents  the  SST  obtained  using  AVHRR 
channels  3  and  4  and  and  7*4  represent  channel  3  and 
channel  4  brightness  temperatures,  respectively. 

A  comparable  empirical  MCSST  algorithm  that  NOAA 
derived  from  drifting  buoy  data  and  used  operationally  for 
NOAA  11  is 

Mr(3,  4)  =  1.586673  “  0.590574 

-E  2.027[sec  (0)  -  1]  +  2.96  (5) 

where  0  is  the  satellite  zenith  angle.  The  term  containing  the 
satellite  zenith  angle  compensates  for  increased  path  length 
for  off-nadir  viewing.  This  correction  for  zenith  angle  allows 
for  comparison  of  satellite  observations  from  several 
AVHRR  swath  angles  with  the  model  results  which,  for 
simplicity,  were  calculated  for  vertical  viewing  only.  The 
close  correspondence  between  channel  coefficients  of  the 
theoretical  algorithm  (4)  and  the  empirical  algorithm  (5) 
indicates  that  the  simple  model  actually  does  a  fair  job  of 
depicting  radiative  transfer  in  the  ocean-atmosphere  system. 

We  now  proceed  to  model  errors  in  MCSST  data.  For  the 
case  of  AVHRR  IR  channels  3  and  4,  we  can  calculate  the 
temperatures  observed  by  the  satellite  using  the  radiative 
transfer  equation  (2).  These  simulated  satellite  observations 
can  then  be  atmospherically  corrected  to  give  7(3,  4)  by 
using  the  theoretical  MCSST  algorithm  (4).  The  resulting 
MCSST  error  E  in  7(3,  4)  can  then  be  calculated; 

£=7(3,  4) -7^  (6) 

or,  after  appropriate  substitution  and  algebraic  manipula¬ 
tion. 


£  =  (1.5t3-0.5t4- 1)(7,-7J  (7) 

Equation  (7)  demonstrates  an  error  in  MCSST  that  is  a 
linear  function  of  the  difference  between  the  sea  surface  and 
atmospheric  temperatures.  7^  is  the  temperature  of  the 
atmosphere  which  has  been  modeled  as  vertically  uniform. 
In  truth,  atmospheric  temperature  varies  with  altitude.  Thus 
(7)  is  difficult  to  compare  with  actual  observations,  as  the 
variable  represents  the  mean  effective  atmospheric  tem¬ 
perature. 

For  each  atmosphere,  LOWTRAN  provides  atmospheric 
temperatures  for  32  levels.  These  temperatures  can  be 
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Fig.  2.  Mean  model  atmospheric  air  temperature  versus  model 
surface  air  temperature  for  five  LOWTRAN  7  atmospheric  models. 
Dotted  line  represents  linear  regression  relationship  obtained  in  (8). 

averaged  over  the  altitudes  representative  of  water  vapor 
attenuation  effects  at  IR  wavelengths  to  obtain  ,  the  mean 
effective  atmospheric  temperature  for  the  spectral  bands 
used.  LOWTRAN  contains  six  atmospheric  models  of  which 
five  are  representative  of  general  atmospheric  conditions 
over  the  ocean  surface  where  SST  retrievjds  are  typically 
obtained  (70°N  to  TO^S).  These  models  are  the  U.S.  Stan¬ 
dard,  tropical,  mid-latitude  summer,  mid-latitude  winter, 
and  Subarctic  summer  atmospheres. 

Atmospheric  model  temperatures  were  averaged  over 
levels  ranging  from  the  surface  to  7  km.  This  altitude  was 
selected  because  (1)  water  vapor  is  the  major  atmospheric 
attenuating  medium  at  the  wavelengths  of  interest  and  (2)  the 
majority  of  atmospheric  water  vapor  typically  exists  in  the 
lowest  7  km,  with  more  than  half  in  the  lowest  3  km  [t/.S.  Air 
Force,  1965]. 

In  Figure  2  the  mean  atmospheric  temperature  is  plotted 
versus  the  model  surface  air  temperature,  which  is  assumed 
to  be  representative  of  the  standard  10-m  height.  Although  a 
limited  number  of  atmospheric  models  are  used,  if  one 
assumes  that  these  models  are  generally  representative  of 
atmospheric  conditions  over  the  global  ocean,  the  results 
demonstrate  a  basic  linear  relationship  between  model  sur¬ 
face  air  temperature  and  mean  model  atmospheric  tempera¬ 
ture.  This  relationship  can  be  expressed  by 

T,,  =  0.874r„(0) 25.86  (8) 

where  7^(0)  represents  the  model  surface  air  temperature. 
The  regressed  correlation  coefficient  between  and  TJO) 
is  excellent  (R  =  0.98).  The  regression  F  statistic  test 
assured  that  the  model  adequately  accounts  for  a  significant 
amount  of  variation  in  the  response  and  that  it  explains  mean 
atmospheric  temperature  at  a  95%  confidence  level. 

Thsi  value  of  from  (8)  can  now  be  substituted  into  (7)  to 
provide  MCSST  errors  in  terms  of  the  sea  surface  tempera¬ 
ture  and  surface  air  temperature. 

£=  (1.5t3 -0.5t4  -  l)[T,  -  0.874r„(0)  -  25.86]  (9) 

Simulation  Experiment 

We  employed  the  LOWTRAN  7  code  to  compute  satellite- 
sensed  apparent  surface  temperatures  for  various  sea  sur¬ 
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face  temperatures  and  atmospheric  profiles  for  the  three 
AVHRR  IR  channels.  The  NOAA  1 1  IR  channel  transfer 
functions  were  used  to  model  the  channel  responses  and 
generate  the  brightness  temperatures  for  each  simulated 
channel.  These  simulated  brightness  temperatures  were  then 
used  to  generate  MCSST,  CPSST,  and  NLSST  retrievals 
under  the  five  atmospheric  model  conditions  mentioned 
previously. 

To  simulate  the  effects  of  ASTD,  the  model  sea  surface 
temperature  was  allowed  to  range  5°C  above  and  below  each 
atmospheric  model  surface  air  temperature  at  increments  of 
rC.  The  error  between  simulated  satellite  SST  algorithm 
retrievals  and  the  sea  surface  temperature  was  then  calcu¬ 
lated  at  each  step,  providing  the  simulated  retrieval  error 
sensitivity  to  ASTD. 

Satellite  SST  retrievals  were  obtained  for  the  three  IR 
channel  combinations  used  operationally;  dual-window 
(channels  3  and  4),  split-window  (channels  4  and  5),  and 
triple-window  (channels  3, 4,  and  5).  The  simulation  retrieval 
algorithms  were  generated  from  multiple  linear  regressions 
of  the  simulated  IR  channel  brightness  temperatures  to 
model  sea  surface  temperatures.  For  simplicity  the  retrieval 
algorithms  were  first  derived  using  cases  in  which  the  sea 
surface  temperatures  equaled  the  model  surface  air  temper¬ 
atures  (ASTD  equal  to  zero).  Algorithms  were  later  derived 
at  ±5°C  ASTD  to  compare  results.  These  results  will  be 
discussed  later.  The  model  MCSST  algorithm  simulation 
equations  for  zero  ASTD  are 

7(3,  4)  =  1.028374+  1.3116(7.,  -  T^)  -  6.89  (10) 

7(3,  4,  5)  =  1.020474  +  0.9464(7,  -  7,)  -  4.73  (11) 

7(4,  5)  =  1.004274  +  3.2014(74  -  7,)  -  0.36  (12) 

The  operational  NOAA  1 1  equation  corresponding  to  7(3. 
4)  was  previously  given  in  (5).  The  following  are  the  opera¬ 
tional  triple-  and  split-window  algorithms  for  NOAA  11. 

MC(3.  4,  5)  =  1.009574  +  1.041(7,  -  7,) 

+  1.76[sec  (0)  -  1]  -  2.04  (13) 

A/C(4.  5)  =  1.056174  +  2.5415(74  -  7,) 

+  0.888(74  -  75)[sec  (0)  -  1]  -  16.9795  (14) 

Differences  between  simulated  and  operational  equations 
may  be  attributed  to  at  least  two  factors:  (1)  the  different 
methods  employed  to  obtain  the  simulated  versus  opera¬ 
tional  equations  and  (2)  the  dependence  of  algorithm  coeffi¬ 
cients  on  ASTD.  The  operational  NOAA  1 1  algorithms  are 
generated  empirically  using  several  hundred  global  drifting 
buoy  SSTs  matched  to  satellite  retrievals  obtained  under 
actual  atmospheric  conditions.  In  comparison,  the  model- 
simulated  algorithms  were  generated  from  a  limited  number 
of  atmospheric  model  temperatures  that  may  or  may  not 
constitute  an  adequate  sample  of  global  conditions.  Since 
water  vapor  continuum  absorption  coefficients  in  the  thermal 
IR  band  are  not  precisely  known,  the  model-simulated 
algorithm  coefficients  strongly  depend  on  the  transmission 
code  used.  The  results  obtained  by  this  study,  however,  will 
show  that  ASTD  error  characteristics  for  model-simulated 
algorithms  and  operational  algorithms  are  very  similar. 

The  other  likely  factor  producing  coefficient  differences 
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will  be  mentioned  here  and  then  expanded  upon  at  the  end  of 
this  section.  It  will  be  demonstrated  that  algorithms  derived 
for  a  specific  mean  ASTD  will  differ  from  algorithms  derived 
at  other  ASTD  conditions  and  will  be  less  accurate  than  the 
other  algorithms  at  those  conditions.  Thus  the  mean  ASTD 
present  in  the  data  set  used  to  derive  the  SST  algorithm  will 
greatly  influence  the  resultant  equation  coefficients  and  the 
equation's  accuracy  at  other  ASTD  conditions. 

Equations  (10),  (1 1),  and  (12)  were  derived  assuming  zero 
ASTD  conditions.  In  reality  the  assumption  that  global  mean 
ASTD  is  zero  may  not  be  true.  This  fact  can  be  easily 
demonstrated  by  utilizing  the  in  situ  SST  and  air  temperature 
data  available  from  the  global  buoy  data  employed  in  this 
study.  This  data  set  reveals  global  buoy  SST  to  be  an 
average  1 .04°C  warmer  than  surface  air  temperature.  Thus 
the  operational  equations  derived  from  global  buoy  data 
would  account  for  whatever  mean  ASTD  exists  in  the  data 
set  used  to  derive  the  equations,  leading  to  a  difference 
between  simulated  and  operational  algorithm  coefficients. 

CPSST  and  NLSST  algorithms  were  also  evaluated.  As 
described  by  Walton  [1988],  the  CPSST  algorithm  assumes  a 
nonlinear  atmospheric  correction  based  upon  IR  brightness 
temperature  magnitude  and  interchannel  temperature  differ¬ 
ence.  The  NLSST  [Walton  et  al.,  1990]  modifies  this  non¬ 
linear  technique  to  be  a  function  of  the  underlying  surface 
temperature.  These  techniques  differ  from  MCSST,  which 
assumes  a  linear  atmospheric  correction  based  solely  upon 
interchannel  temperature  differences.  The  simulated  and 
operational  NOAA  11  CPSST  and  NLSST  algorithms  are 
listed  in  the  appendix.  The  T f  term  in  the  NLSST  equations 
represents  the  underlying  surface  temperature  obtained  from 
an  analysis  field,  climatology,  or  MCSST  retrieval.  This 
study  used  the  latter  value  for  T/ .  applying  the  appropriate 
MCSST  algorithm  (5,  13,  or  14). 

Figure  3  depicts  the  simulated  MCSST  retrieval  error  for 
the  three  MCSST  retrieval  algorithms  in  each  of  the  model 
atmospheres.  The  dual-,  split-,  and  triple-window  algorithm 
errors  are  denoted  by  solid,  dotted,  and  dashed  lines, 
respectively.  Each  algorithm  was  derived  assuming  zero 
ASTD,  and  it  is  interesting  to  note  that  with  few  exceptions, 
each  is  most  accurate  under  such  conditions.  Also,  note  that 
the  simulated  dual-window  and  triple-window  algorithms 
consistently  overestimate  SST  when  the  air  temperature  is 
much  warmer  than  the  underlying  sea  surface  and  underes¬ 
timate  SST  when  conditions  are  reversed.  The  split-window 
algorithm  shows  much  less  consistent  sensitivity  to  ASTD. 

The  slopes  for  these  error  lines  are  presented  in  Table  1 . 
These  values  were  obtained  via  a  linear  least  squares  fit  of 
the  error  at  each  TC  increment.  Slope  t  statistics  were 
calculated  for  each  fit  at  the  95%  confidence  level  to  test  the 
hypothesis  that  the  true  slope  may  actually  equal  zero. 
Those  slope  values  with  t  statistics  not  found  to  be  statisti¬ 
cally  significant  are  enclosed  in  parentheses.  The  hypothesis 
that  the  slope  is  actually  zero  cannot  be  rejected  for  such 
cases.  Note  that  only  two  model  split-window  slopes  can  be 
identified  statistically  as  nonzero.  The  overall  trend  for  the 
other  two  algorithms  is  a  statistically  significant  slope  that 
varies  by  only  a  small  amount  between  atmospheric  models. 

Both  the  dual-window  and  triple-window  algorithms  indi¬ 
cate  significant  MCSST  retrieval  error  in  all  model  atmo¬ 
spheres  as  the  difference  between  the  sea  surface  tempera¬ 
ture  and  the  air  temperature  incrcu.3e$.  This  finding  predicts 
the  MCSST  retrieval  to  be  colder  than  the  true  sea  surface 


temperature  when  the  air  temperature  is  significantly  less 
than  the  sea  surface  temperature.  Similarly,  the  MCSST 
retrieval  is  predicted  to  be  warmer  than  the  true  sea  surface 
temperature  when  the  air  temperature  is  substantially 
greater  than  the  sea  surface  temperature.  Taking  an  average 
of  the  simulated  slope  values,  the  dual-window  algorithm 
error  is  approximately  1°C  for  every  7°C  change  in  the 
ASTD.  For  the  triple-window  algorithm  the  simulation  pre¬ 
dicts  a  1°C  error  for  every  10°C  change  in  the  ASTD. 

The  split-window  algorithm  only  shows  significant  slope 
for  tropical  and  mid-latitude  winter  atmospheres  with  no 
statistically  significant  slope  for  the  other  atmospheric  mod¬ 
els.  These  predictions  indicate  that  the  dual-window  and 
triple-window  MCSST  algorithms  should  be  more  sensitive 
to  the  ASTD  than  is  the  split-window  algorithm.  This  finding 
is  generally  expected,  since  the  MCSST  algorithm  assumes 
that  the  atmospheric  absorption  process  is  similar  at  the  two 
wavelengths  used  for  the  atmospheric  correction.  Such 
assumptions  are  obviously  more  valid  for  channels  4  and  5 
[Prabhakara  et  al.,  1974],  which  are  centered  within  1  ton  of 
each  other. 

Tables  2  and  3  display  the  predicted  slopes  for  CPSST  and 
NLSST  errors,  respectively.  Although  the  dual-  and  triple¬ 
window  error  slopes  are  comparable  to  the  respective  MC¬ 
SST  algorithm  errors,  both  the  CPSST  and  NLSST  split- 
window  algorithms  demonstrate  greater  sensitivity  to  the 
ASTD  than  does  the  MCSST  split-window  algorithm.  These 
split-window  algorithms  may  be  in  error  TC  for  every 
I0®-I2°C  change  in  ASTD.  These  results  differ  from  the 
MCSST  split-window  results  which  do  not  show  consistent 
sensitivity.  This  difference  may  be  due  to  the  increased 
dependence  of  the  CPSST  algorithm  on  the  brightness  tem¬ 
perature  magnitude  and  to  the  NLSST  algorithm  depen¬ 
dence  on  the  underlying  surface  temperature  magnitude.  Not 
surprisingly,  this  dependence  makes  the  algorithms  more 
sensitive  to  ASTD  changes. 

As  stated  previously,  simulation  equations  were  also 
derived  under  model-simulated  conditions  in  which  the 
ASTD  was  nonzero,  specifically,  ±5^.  The  resulting  algo¬ 
rithms  generated  error  slopes  nearly  identical  to  the  respec¬ 
tive  algorithms  in  Tables  1-3.  However,  the  error  magni¬ 
tudes,  or  bias,  of  each  algorithm  were  displaced 
approximately  0.5°C  on  either  side  of  the  slope  lines  gener¬ 
ated  from  the  zero  ASTD  algorithms. 

An  example  of  this  is  shown  in  Figure  4  for  the  U.S. 
Standard  Atmosphere  model  dual-window  algorithm.  This 
figure  clearly  shows  that  an  algorithm  derived  under  one 
mean  ocean-atmosphere  relationship  will  not  be  as  accurate 
under  other  ocean-atmosphere  conditions.  Each,  however, 
is  most  accurate  when  ASTD  equals  the  mean  conditions 
used  to  derive  it.  Other  simulated  retrieval  algorithms  and 
model  atmospheres  generated  similar  results  regarding  slope 
and  algorithm  error  bias.  This  difference  demonstrates  the 
importance  of  the  mean  air-sea  temperature  contrast  used  to 
derive  each  simulated  algorithm. 

This  result  agrees  with  the  findings  of  Eyre  [1987]  which 
demonstrated  that  retrieval  errors  are  minimized  for  physical 
ocean-atmosphere  situations  that  equal  the  mean  ocean- 
atmosphere  state  used  to  derive  the  algorithm.  An  algorithm 
derived  under  zero  ASTD  will  perform  best  for  such  condi¬ 
tions.  Likewise,  an  algorithm  derived  from  ASTD  of  5“C  will 
be  less  accurate  than  the  first  algorithm  at  zero  ASTD  but 
more  accurate  at  5°C  difference. 
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Fig.  3.  Simulated  MCSST  retrieval  error  as  a  function  of  SST-air  temperature  difference  by  model  atmosphere  for 
each  of  three  nighttime  algorithms.  Dual  window  (channels  3  and  4)  is  denoted  by  a  solid  line.  Split  window  (channels 
4  and  5)  is  denoted  by  a  dotted  line.  Triple  window  (channels  3, 4,  and  5)  is  denoted  by  a  dashed  line,  (a)  Tropical,  {b) 
mid-latitude  summer,  (c)  mid-latitude  winter,  (<f)  Subarctic  summer,  and  (e)  U.S.  Standard. 


Although  algorithms  can  be  generated  to  reduce  the  error 
bias  at  specific  mean  ocean-atmosphere  states,  the  error  of 
each  algorithm  is  still  sensitive  to  changes  in  the  actual 
ASTD.  This  is  clearly  evident  in  Figure  4.  Each  algorithm  is 
most  accurate  at  the  mean  ASTD  used  to  derive  it.  How¬ 
ever,  each  algorithm  error  also  follows  approximately  the 
same  slope  with  respect  to  ASTD  changes.  Thus  ASTD 
changes  will  have  similar  effects  on  regionally  as  well  as 


globally  derived  algorithm  accuracies.  Such  effects  will  be 
negligible  only  if  ASTD  varies  less  than  a  few  degrees  from 
mean  conditions. 

Ground  Truth  Comparisons 

Operational  satellite  retrievals  are  matched  in  near-real 
time  to  colocated  buoy  measurements  as  part  of  the  NESDIS 
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TABLE  1.  Predicted  MCSST  Error  Slopes  TABLE  3.  Predicted  NLSST  Error  Slopes 


Atmospheric  Model 

TO.  4) 
Slope 

7-(4.  5) 
Slope 

T(3,  4.  5) 
Slope 

Atmospheric  Model 

r(3,  4) 
Slope 

r(4. 5) 
Slope 

r(3,  4,  5) 
Slope 

Tropical 

-0.184 

0.079 

-0.103 

Tropical 

-0.142 

-0.045 

-0.070 

Midlatitude  summer 

-0.168 

(0.010) 

-0.114 

Midlatitude  summer 

-0.166 

-0.129 

-0.130 

Midlatitude  winter 

-0.097 

-0.041 

-0,081 

Midlatitude  winter 

-0.107 

-0.128 

-0.090 

Subarctic  summer 

-0.136 

(0.020) 

-0.090 

Subarctic  summer 

-0.148 

-0.125 

-0.120 

U.S.  Standard 

-0.140 

(0.009) 

-0.095 

U.S.  Standard 

-0.132 

-0.112 

-0.107 

Average 

-0.145 

0.015 

-0.097 

Average 

-0.139 

-0.108 

-0.103 

operational  retrieval  generation  process.  These  matches  are 
retained  in  a  matchup  data  base  that  can  be  used  for 
subsequent  off-line  analysis.  Colocated  satellite  and  buoy 
SST  matchups  were  obtained  from  this  data  base  for  the 
months  of  February  through  May  1991.  Matchups  with 
moored  and  drifting  buoys  were  collected  globally  to  provide 
a  diverse  sample  of  matchups  in  all  seasons. 

For  an  ideal  comparison  the  in  situ  and  spacebome 
sensors  should  measure  the  same  quantity,  which  would 
require  the  use  of  a  near- surface  infrared  radiometer.  Since 
buoys  typically  measure  SST  at  a  depth  of  1  m,  it  is  argued 
that  these  “bulk  SST”  measurements  can  be  decoupled  from 
the  true  surface  temperature  by  near-surface  temperature 
gradients,  the  so-called  “skin  effect”  [Robinson  et  al.,  1984; 
Hepplewhite,  1989J. 

Given  that  infrared  radiometer  in  situ  measurements  from 
research  ships  are  quite  scarce  and  very  regionalized  [e.g., 
Schluessel  et  al.,  1987;  Coppin  et  al.,  1991]  and  that  atmo¬ 
spheric  and  oceanographic  models  and  analyses  traditionally 
assimilate  bulk  SST,  it  is  very  convenient  to  regress  the 
satellite  retrievals  to  this  more  readily  available  parameter. 
Use  of  a  regionally  derived  in  situ  radiometer  algorithm  as  a 
global  algorithm  can  also  be  questioned  because  of  the 
limited  spatial  and  seasonal  extent  of  the  matchup  compar¬ 
isons  used  to  generate  the  algorithm.  Buoys  provide  a 
globally  diverse  set  of  observations  accurate  to  O.TC  that 
can  be  readily  used  to  generate  global  algorithms  and 
monthly  global  accuracy  statistics. 

In  order  to  limit  the  time  and  distance  variability  that  can 
occur  between  buoy  and  satellite  data  the  matchups  used  in 
this  study  include  only  those  buoy  observations  within  1 
hour  and  10  km  of  each  satellite  retrieval.  Minnett  [1991] 
demonstrates  that  such  time  and  distance  criteria  provide 
representative  matchups  by  eliminating  possible  temporal 
and  spatial  differences  that  might  occur  between  satellite  and 
in  situ  data  because  of  diurnal  SST  changes  or  ocean  frontal 
gradients. 

Because  some  satellite  retrievals  are  matched  to  more 
than  one  buoy  observation,  the  data  set  was  rehned  to 
include  just  the  single  closest  buoy  matchup  in  time  and 


TABLE  2.  Predicted  CPSST  Error  Slopes 


Atmospheric  Model 

7'(3,  4) 
Slope 

r(4, 5) 
Slope 

r(3,  4.  5) 
Slope 

Tropical 

-0.155 

-0.053 

-0.087 

Midlatitude  summer 

-0.158 

-0.090 

-0.106 

Midlatitude  winter 

-0.1.32 

-0.044 

-0.060 

Subarctic  summer 

-0.145 

-0.058 

-0.082 

U.S.  Standard 

-0.133 

-0.056 

-0.080 

Average 

-0.145 

-0.060 

-0.083 

distance  to  each  retrieval.  The  matchups  were  further  re¬ 
stricted  to  nighttime-only  satellite  retrievals  to  avoid  diurnal 
effects.  The  resulting  data  set  includes  160  matchups  ob¬ 
tained  from  most  ocean  regions  of  the  world  (Figure  5). 

Using  (5),  AVHRR  channels  3  and  4  brightness  tempera¬ 
tures  were  used  to  calculate  MC  (3,  4)  for  each  of  the  160 
satellite-buoy  matchups.  Figure  6  shows  the  NESDIS  dual¬ 
window  MCSST  retrievals  minus  buoy  SST  versus  the  buoy 
SST  minus  air  temperature.  Note  the  negative  slope  of  a 
satellite  minus  buoy  values  evident  in  Figure  6  u*  the  buoy 
SST  minus  air  temperature  values  increase.  A  simple  least 
squares  regression  fit  to  the  data  results  in  a  satellite  minus 
buoy  error  slope  of  -0. 143  with  a  standard  error  of  0.03  at 
the  95%  confidence  level. 

To  test  whether  the  true  value  of  the  slope  is  equal  to  zero, 
the  critical  t  statistic  is  calculated  to  be  1 .645  for  158  degrees 
of  freedom  at  the  95%  confidence  level.  The  regression 
model  slope  t  statistic  is  found  to  be  -3.507,  which  rejects 
the  hypothesis  that  the  true  slope  is  equal  to  zero.  The 
regression  slope  is  therefore  found  to  be  statistically  signif¬ 
icant.  The  average  slope  of  -0. 145  predicted  from  the  model 
atmospheres  is  within  the  standard  error  of  the  in  situ  data 
slope  estimate.  These  encouraging  results  confirm  the  simu¬ 
lation  experiment  predictions  of  dual-window  retrieval  er¬ 
rors  approaching  TC  for  every  7°C  change  in  the  ASTD. 

Figure  7  displays  the  NESDIS  triple-window  MCSST 
algorithm  retrievals  minus  buoy  SST  versus  the  buoy  SST 
minus  air  temperature.  These  results  were  obtained  using 


SST  -  AIR  TEMPERATURE  fC) 

Fig.  4.  Simulated  dual-window  MCSST  retrieval  error  as  a 
function  of  SST-air  temperature  difference  for  the  U.S.  Standard 
Atmosphere  model.  The  three  lines  represent  errors  associated  with 
model-simulated  MCSST  algorithms  derived  under  three  different 
air-sea  temperature  conditions:  plus  5  (+5).  zero  (0),  and  minus  5 
(-5). 
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Fig.  5.  Location  of  buoy  in  situ  data  used  in  data  analysis.  Each  satellite-buoy  matchup  location  is  denoted  by  a  plus. 


(13).  Like  the  dual-window  algorithm,  a  negative  slope  is 
evident  in  the  figure.  A  simple  least  squares  regression  fit  to 
the  data  results  in  a  satellite  minus  buoy  error  slope  of 
-0.101  with  a  standard  error  of  0.03.  The  regression  model 
slope  t  statistic  is  calculated  to  be  -3.426,  which  is  statisti¬ 
cally  significant  at  the  95%  confidence  level.  This  slope 
agrees  quite  well  with  the  average  slope  of  -0.097  predicted 
by  the  simulation  experiment.  These  results  predict  satellite 
retrieval  errors  approaching  TC  for  ASTD  changes  of  10°C. 

Figure  8  shows  the  NESDIS  split-window  MCSST  algo¬ 
rithm  retrievals  minus  buoy  SST  versus  the  buoy  ASTD. 
These  results  were  obtained  using  (14).  This  figure  indicates 
a  more  random  distribution  of  matchups  with  ASTD  than 
does  the  previous  figures.  A  simple  least  squares  regression 
fit  predicts  a  slope  of  -0.043  with  a  standard  error  of  0.044. 
The  slope  t  statistic  is  calculated  as  -0.%9,  which  is  not 
statistically  significant  at  the  95%  confidence  level.  Thus  the 


BUOY  SST  -  AIR  TEMPERATURE  (‘C) 


Fig.  6.  MCSST  dual-window  retrieval  minus  buoy  SST  plotted 
versus  the  buoy  SST  minus  buoy  air  temperature. 


hypothesis  that  the  true  slope  value  is  zero  cannot  be 
rejected.  This  result  is  in  general  agreement  with  the  simu¬ 
lation  experiment,  which  indicated  no  significant  slope  with 
ASTD  except  for  the  tropical  atmosphere. 

Table  4  lists  the  error  slopes  for  each  algorithm  type  when 
compared  to  the  buoy  matchup  data  set.  Similar  to  the 
MCSST  algorithms,  the  average  CPSST  and  NLSST  versus 
buoy  error  slopes  were  found  to  be  within  the  standard  error 
of  each  slope  predicted  by  the  model  atmospheres.  These 
matchup  results  agree  favorably  with  the  simulated  predic¬ 
tions,  showing  the  dual-window  and  triple-window  algo¬ 
rithms  to  be  more  sensitive  to  ASTD  than  is  the  split-window 
algorithm,  especially  for  the  MCSST  technique. 

Table  4  also  shows  that  the  nonlinear  CPSST  and  NLSST 
split-window  algorithms  are  more  sensitive  to  ASTD  than  is 
the  MCSST  split- window  algorithm.  This  difference  was  also 
demonstrated  in  the  model  predictions.  Thus  both  model 


Fig.  7.  MCSST  triple-window  retrieval  minus  buoy  SST  plotted 
versus  the  buoy  SST  minus  buoy  air  temperature. 
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BUOY  SST  -  AIR  TEMPERATURE  ('C) 

Fig.  8.  MCSST  split-window  retrieval  minus  buoy  SST  plotted 
versus  the  buoy  SST  minus  buoy  air  temperature. 


predictions  and  ground  truth  comparisons  indicate  that  use 
of  satellite  SST  retrievals  under  extremely  anomalous  con¬ 
ditions  of  ASTD  can  generate  significant  retrieval  errors. 

Table  S  lists  the  mean  bias  and  root  mean  square  differ¬ 
ence  (rmsd)  for  each  operational  algorithm  as  compared  to 
global  buoy  SSTs.  The  table  displays  global  as  well  as 
regional  statistics.  The  statistics  have  also  been  divided  into 
atmospheric  moisture  classes,  using  the  channel  4  minus  S 
values  as  a  measure  of  atmospheric  moisture  content  and  the 
buoy  SST  minus  air  temperature  (BSST  minus  AT)  to 
classify  statistics  by  air-sea  interface  conditions. 

In  general,  large  positive  BSST  minus  AT  corresponds  to 
large  negative  retrieval  bias.  Similarly,  negative  or  very  low 
positive  BSST  minus  AT  corresponds  to  positive  or  very  low 
negative  retrieval  bias.  This  supports  the  model  results 
obtained  earlier,  demonstrating  that  when  air  temperature  is 
significantly  greater  (less)  than  SST,  the  satellite  retrieval 
will  overestimate  (underestimate)  SST 

This  result  is  not  evident  in  all  cases,  particularly  the 
split-window  MCSST  retrievals  which  are  less  sensitive  to 
ASTD.  The  biases  observed  in  the  MCSST  split-window 
algorithm  regional  cells  are  more  a  result  of  sensitivity  to  the 
4  minus  5  values  or  atmospheric  moisture  content.  Some  of 
the  largest  negative  biases  in  Table  5  are  present  in  these 
cells.  It  appears  that  the  split-window  retrievals  are  too  cold 
when  the  channel  4  minus  5  values  are  very  small  and  too 
warm  when  conditions  are  reversed. 

The  nonlinear  CPSST  and  NLSST  algorithms  improve  on 
the  MCSST  split-window  algorithm  error  quite  nicely,  re¬ 
ducing  bias  and  rmsd  in  almost  every  region.  These  statistics 
demonstrate  the  improved  retrieval  accuracy  that  can  be 
obtained  using  a  nonlinear  split-window  algorithm.  The 
exception  to  this  is  when  the  air  temperature  exceeds  SST 


TABLE  4.  Buoy  Matchup  Error  Slopes 


SST 

Algorithm 

T(3,  4) 

Slope 

T(4,  5) 

Slope 

r(3,  4,  5) 
Slope 

MCSST 

-0.143 

(-0.043) 

-O.IOt 

CPSST 

-0.136 

-0.84 

-0.105 

NLSST 

-0.143 

-0.106 

-0.129 

and  when  SST  substantially  exceeds  air  temperature.  These 
results  demonstrate  the  increased  sensitivity  of  the  CPSST 
and  NLSST  split- window  algorithms  to  ASTD. 

Comparison  between  the  dual-,  split-,  and  triple-window 
techniques  shows  that  the  triple-window  algorithms  demon¬ 
strate  better  overall  rmsd  statistics  than  the  others.  The 
triple-window  algorithm  bias  statistic,  however,  is  signifi¬ 
cantly  sensitive  to  ASTD  changes.  The  dual-window  algo¬ 
rithms  also  demonstrate  significant  bias  sensitivity  as  ASTD 
changes.  This  sensitivity  is  somewhat  greater  than  the  triple¬ 
window  algorithm,  as  evidenced  in  the  rows  depicting  ASTD 
conditions. 

Overall,  none  of  the  160  BSST  minus  AT  matchups 
exceeded  5°C  or  were  less  than  -2°C.  This  limited  range  of 
difference  is  related  to  the  general  world  o''ean  situation  in 
which  the  majority  of  ASTDs  are  less  than  ±2°C.  The 
formation  of  cloud  or  fog  as  the  temperature  difference 
reaches  one  extreme  or  the  other  may  also  impair  limits  on 
these  upper  and  lower  thresholds.  Thus,  although  such 
conditions  could  potentially  cause  significantly  erroneous 
satellite  SST  retrievals,  cloud  formation  would  prevent  the 
generation  of  satellite  retrievals. 

Derivation  of  regionally  exclusive  retrieval  algorithms 
may  not  be  an  adequate  solution  to  the  ASTD  sensitivity  of 
retrievals.  ASTD  is  a  continually  changing  quantity,  and  one 
regional  algorithm  will  not  always  serve  a  specific  region's 
range  of  ASTD  conditions  unless  ASTD  never  varies  more 
than  a  few  degrees.  Use  of  air  temperature  information  could 
possibly  provide  improved  retrieval  accuracy  when  large 
ASTD  does  exist.  Air  temperature  information  would  need 
to  have  a  fine  enough  time  and  spatial  resolution  that  its  use 
would  not  generate  errors  greater  than  those  being  cor¬ 
rected.  Operational  applications  would  be  best  served  by  air 
temperature  information  that  is  both  temporally  (within  1 
hour)  and  spatially  (within  10-20  km)  accurate.  However, 
given  the  present  spatial  resolution  of  global  analyses  and 
the  relatively  rapid  changes  that  can  occur  in  air  tempera¬ 
ture,  operational  use  of  such  information  is  conjectured  to  be 
inadequate  for  present  needs. 

Discussion  and  Conclusions 

A  simple  model  has  been  developed  that  predicts  the 
dependence  of  MCSST  satellite  retrieval  errors  on  the 
ASTD.  A  simulation  experiment  has  been  performed  to  test 
this  dependence.  Three  satellite  retrieval  algorithms,  the 
dual-window  (channels  3  and  4),  the  split-window  (channels 
4  and  5),  and  the  triple-window  (channels  3,  4,  and  5),  plus 
the  nonlinear  CPSST  and  NLSST,  were  simulated  and 
compared.  Simulation  results  have  also  been  compared  to 
global  buoy  SST  and  air  temperature  data  matched  to 
satellite  retrievals  within  specific  time  and  distance  con¬ 
straints.  It  is  evident  from  the  comparisons  that  the  model 
results  simulate  the  effect  of  ASTD  present  in  the  satellite- 
buoy  matchups  and  that  the  use  of  satellite  SST  retrievals 
under  anomalous  ASTD  conditions  can  generate  significant 
retrieval  errors. 

In  general,  when  the  actual  ASTD  is  close  to  the  mean 
ASTD  conditions  present  in  the  algorithn.  derivation  data 
set,  the  satellite  retrieval  algorithm  will  approximate  ground 
truth  SST  measurements  quite  well.  However,  it  is  also 
revealed  that  global  satellite  retrievals  can  overestimate  SST 
when  the  SST  minus  air  temperature  difference  significantly 


TABLE  S.  Sateliile/Buoy  Matchup  Statistics  for  the  World  and  Various  Global  Regions 
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exceeds  the  mean  difference  present  in  the  algorithm  deri¬ 
vation  data  set.  Likewise,  retrievals  can  underestimate  SST 
when  conditions  are  reversed.  Dual-window  satellite  retriev¬ 
als  may  be  in  error  1°C  for  every  7°C  change  in  ASTD  from 
the  mean  conditions.  Likewise,  triple-window  retrievals  may 
be  in  error  TC  for  every  10“C  change  in  ASTD. 

The  split-window  retrievals  demonstrate  less  consistent 
sensitivity  to  ASTD  except  for  the  nonlinear  CPSST  and 
NLSST  algorithms.  The  lack  of  sensitivity  to  ASTD  in  the 
MCSST  method  is  most  likely  a  result  of  the  linear  atmo¬ 
spheric  correction  technique  and  the  proximity  of  channels  4 
and  5  (10.5  and  1 1.5  ^m).  Since  channel  3  is  centered  at  3.7 
lun,  the  assumption  that  atmospheric  absorption  is  similar  in 
each  IR  channel  is  less  valid  for  the  dual-  and  triple-window 
algorithms  than  for  the  split- window  algorithm.  The  nonlin¬ 
ear  CPSST  and  NLSST  split-window  algorithms  demon¬ 
strate  greater  dependence  on  the  channel  brightness  temper¬ 
ature  and  underlying  surface  temperature  magnitudes, 
resulting  in  greater  sensitivity  to  ASTD  changes.  Nonlinear 
split- window  retrieval  algorithm  errors  of  TC  can  occur  for 
10°-12°C  changes  in  ASTD. 

In  conclusion,  satellite-derived  SSTs  in  ocean  regions 
susceptible  to  extreme  ASTD  conditions  may  be  suspect 
when  using  an  algorithm  derived  under  less  extreme  condi¬ 
tions.  Such  circumstances  could  generate  significant  satellite 
retrieval  errors.  Equations  can  be  derived  to  attain  better 
accuracy  for  extreme  atmospheric  situations;  however, 
these  equations  would  demonstrate  significant  error  when 
air-sea  temperature  conditions  become  less  extreme.  Thus 
selection  of  regionally  exclusive  algorithms  may  not  be  an 
adequate  solution,  because  ASTD  is  nonstatic,  and  one 
regional  algorithm  will  not  always  serve  that  specific  region’s 
range  of  ASTD  conditions  unless  the  range  is  less  than  a  few 
degrees. 

Retrieval  accuracy  could  theoretically  improve  by  incor¬ 
porating  regionally  accurate  air  temperature  information  into 
the  algorithm.  Operational  use  of  currently  available  global 
analysis  information  is,  however,  probably  inadequate  ow¬ 
ing  to  analysis  spatial  resolution  and  the  relatively  rapid  air 
temperature  changes  that  can  occur  between  analysis  prod¬ 
ucts.  Satellite  SST  retrieval  users  should  thus  be  aware  of 
both  the  mean  atmosphere-ocean  state  used  to  derive  the 
SST  algorithm  in  use  and  the  ASTD  conditions  to  which  it  is 
applied.  Results  demonstrate  that  satellite  retrieval  errors  of 
O.r  to  0.14°C  are  possible  for  every  degree  that  ASTD 
differs  from  the  mean  atmosphere-ocean  state  used  to  derive 
the  algorithm. 


Appendix:  Simulated  and  Operational  CPSST 
AND  NLSST  Algorithms 


Simulated 


CP(l,  4)  =  0.995 IT4  -I- 


O.O736T4  -  19.5579 
0.07414  -  0.0383r3  -  9.0576 


0.097 ITj  -  26.04 
.09757,  -  0.07474  -  6.26 

•  (74  -Ts  +  0.239)  -  271.89 
NL(3,  4)  =  1.031173  -I-  0.009257/(73  "  T*)  “  280.84 
NL(3,  4,  5)  =  1.023874  +  0.7351(73  -  7,) 

-F  0.006827/(73  -  7,)  -  278.87 
NL(4,  5)  =  1.030274  -F  0.8884  7/(74  -  7,)  -  280.73 
Operational 


CP(4,  5)  =  0.99557,  -F 


CP(3,  4)  =  0.985374  + 


0.170874  -  58.47 


\^0. 1733474  -  0.07757,  -  33.74 

•(7,  -  74  -  6.44)  -F  1.97[sec  (»)-!]-  257.28 

CP(3,  4,  5)  =  0.971274 

/  0.168474  -  34.32  \ 

^  \0.2C527,  -  0.07757,  -  20.01  ~ 

+  1.87(sec  (e)  -  1)  -  276.59 

0.19607,  -  48.61 


CP(4,  5)  =  0.95487,  + 


\0.20527, -0.173374-6.11^ 

•  (74  -  7,  -F  1.46)  -F  0.98(74  -  7,) 
•,sec(e)-  1]-  263.84 
NL(3,  4)  =  1.0087,  -F  0.021167/(7,  -  74) 

-F  2.065[sec  (0)  -  1]  -  273.69 
NL(3,  4,  5)  =  1.000674  -F  0.245(7,  -  7,) 

-F  0.027667/(7,  -  7,)  +  1.88[sec  (»)  -  1]  -  272.36 
NL(4,  5)  =  0.960474  +  0.087527/(74  -  7,) 

-F  0.852(74  -  7,)[sec  (6)  -  1]  -  261.46 
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